On the theory of resonant susceptibility of dielectric glasses in magnetic field 
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The anomalous magnetic field dependence of dielectric properties of insulating glasses in the tem- 
perature interval WmK < T < 50mA' is considered. In this temperature range, the dielectric 
permittivity is defined by the resonant contribution of tunneling systems. The external magnetic 
field regulates nuclear spins of tunneling atoms. This regulation suppresses a nuclear quadrupole in- 
teraction of these spins with lattice and, thus, affects the dielectric response of tunneling systems. It 
is demonstrated that in the absence of an external magnetic field the nuclear quadrupole interaction 
b results in the correction to the permittivity 8\ ~ \b\ /T in the temperature range of interest. An 
application of a magnetic field results in a sharp increase of this correction approximately by a factor 
of two when the Zeeman splitting m approaches the order of |6|. Further increase of the magnetic 
field results in a relatively smooth decrease in the correction until the Zeeman splitting approaches 
the temperature. This smooth dependence results from tunneling accompanied by a change of the 
nuclear spin projection. As the magnetic field surpasses the temperature, the correction vanishes. 
The results obtained in this paper are compared with experiment. A new mechanism of the low 
temperature nuclear spin-lattice relaxation in glasses is considered. 

PACS numbers: 6143.Fs, 77.22.Ch,75.50Lk 
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I. INTRODUCTION 

In 1998, Strehlow et al. 1 discovered and anomalous 
low-temperature sensitivity of the dielectric properties 
of various dielectric glasses to a magnetic field. This 
findings continue to attract the attention both of exper- 
imentalists and theoristsi 2 i 3 i 4 i 5 i 6 i 7 i 8 i 9 i 10 i 1:L As a result of 
further investigations, it was found that the effect was 
anomalously pronounced for an insulator in the absence 
of magnetic impurities. 

Many low-temperature properties of glasses are suc- 
cessfully described by a standard tunneling model. 12 ::iJ ! 
In a good approximation the tunneling systems (TS) can 
be treated as a particle moving in a double- well potential 
(DWP). After observing of the anomalous glassy behav- 
ior in a magnetic field, several extensions of the standard 
tunneling model have been suggested— In our opin- 
ion, the model proposed by A. Wurger, A. Fleischmann, 
C. Enssii is the most viable. It assumes that the tunnel- 
ing particle possesses a nuclear electric quadrupole mo- 
ment. As a result, the particle energy acquires an extra 
splitting b in the crystal electric field gradient (EFG). It 
is important that, in general case, the local axis of EFG 
differ in different wells of DWP (see. Fig. The mag- 
netic field then interacts with the nuclear spin magnetic 
moment and results in the Zeeman splitting. This mod- 
ifies nuclear spin states in each well thus affecting the 
tunneling system properties. 

The echo experiments provide convincing evidence of 
the influence of the nuclear quadrupole moments on 
tunnclingiSiSiS Recently Nagel et al— investigated an 
isotope effect in polarization echo experiments in glasses. 
They observed the qualitatively different echo spectra on 
amorphous glycerol when hydrogen, which has a zero 



electric quadrupole moment, was substituted by deu- 
terium, which possesses a nonzero quadrupole moment. 
Thereupon, one can conclude that the quadrupole elec- 
trical moment of the TS's is the key feature responsible 
for the effect. 

To our knowledge, there are several theoretical stud- 
ies of the nuclear quadrupole effect on glassy dielec- 
tric propertiesiiSiiSiSSi^iiS 2 ! In papers22i2i we have inves- 
tigated the low temperature limit T < b. We have shown 
that in this regime the nuclear quadrupole interaction re- 
sults in the breakdown of coherent tunneling due to the 
formation of a nuclear spin polaron state of tunneling 
particle. It was demonstrated that for T < lOmK the 
saturation in the dielectric constant temperature depen- 
dence observed in a variety of glasses can be explained 
by that breakdown. 

The case of high temperatures T ^> b has been ex- 
amined in Refs— *i2i— Paper— is concerned with the 
many - body relaxation of tunneling system experienc- 
ing quadrupole and Zeeman splittings. Ref— addresses 
a nonlinear behavior of dielectric permittivity at high 
temperatures. Paper— investigates the correction to 
the real part of the permittivity Sx due to the nuclear 
quadrupole interaction and the magnetic field. The au- 
thors of paper— have supposed that this correction is due 
to tunneling systems, which have energy splitting of the 
order of temperature. As a result, a very weak contribu- 
tion Sx has been found to possess a strong temperature 
dependence Sx ~ 1/T 6 , which differs significantly from 
the experimental behavior of Sx ~ 

In the present paper we argue that in the case \b\ <C 
T, only tunneling systems are important, which have 
energy splitting of the order of quadrupole splitting 
A, Ao ~ b <C T. It is shown that for vanishing magnetic 
field the correction to the standard dielectric permittivity 
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X can be expressed as Sx/x ~ H After application of 
magnetic field, Sx first increases sharply. When the Zee- 
man splitting passes the quadrupole splitting, Sx reaches 
the maximum followed by a relatively smooth decrease, 
yet, the dependence Sx ~ \b\ /T remains. The suggested 
theory agrees qualitatively with the experimental data^S 
The paper is organized as follows. In Section [H] we 
introduce the tunneling model with the quadrupole and 
Zeeman interactions. As an example, we derive the tun- 
neling Hamiltonian for the nuclear spin 7 = 1. Then, in 
Sec. IIIII we investigate the correction to the resonant 
permittivity induced by the quadrupole and the Zeeman 
interactions. In the following sections we present numer- 
ical analysis for the resonant permittivity in the case of 
nuclear spin 1=1, For the cases 1 = 1 and I = 3/2, we 
propose an exact analytical solution for the permittivity 
Sx induced only by the quadrupole interaction. We give 
the qualitative explanation for the temperature and the 
magnetic field dependencies. In conclusion, we consider 
the relation of the results obtained to the experimental 
data. Also, a new mechanism of the low temperature 
nuclear spin-lattice relaxation in glasses is considered, 
which is due to tunneling system quadrupole interaction. 



II. GENERALIZATION OF THE STANDARD 
TUNNELING MODEL 

According to the standard tunneling modeii2ii&ii at 
low-temperatures, amorphous solids are represented by 
an ensemble of tunneling systems. One can conceive a 
tunneling system as a particle moving in a double-well 
potential (DWP) characterized by the asymmetry energy 
A and the tunneling splitting Aq . Motion of such a parti- 
cle is described by the standard two-level pseudospin 1/2 
Hamiltonian 



n = ——o —a 

2 2 



(1) 



Following Refs.jiSiiii^ we assume that parameters A, Ao 
obey the universal distribution 



P(A,A„) = ^ 



(2) 



where P is a constant. 

The tunneling particle possesses its own internal de- 
grees of freedom associated with its nuclear spin I. The 
tunneling particle gains Zeeman energy in the magnetic 
field B 



E int = -g/3BI, 



(3) 



where g is the Lande factor and (3 is the nuclear magne- 
ton. Typically the product g(3B reaches the value lmK 
at B w 5T. 

Consider the case of the spin I > 1. In this case the 
nucleus can possess an electric quadrupole moment. It 
interacts with the crystal field characterized by the tensor 




FIG. 1: Tunneling system with asymmetry energy A, and 
the corresponding quadrupole quantization axes in the left 
well ul and in the right well ur (see Ref^i) 



of electric field gradient (EFG) qij. The Hamiltonian 
of the particle interacting with the crystal field can be 
expressed as 2 ^ 



Hr 



b if 



I 2 /3 



(4) 



Here, the parameter b 



3e Qgn 



41(21-1) designates the 
quadrupole interaction constant and the asymmetry pa- 
rameter is given by 



g22 — 933 
9n 



(5) 



We assume that the Cartesian ei,e2,es axes are chosen 
so that < 922 < 9n, since then < je < 1. If h = 0, 
then the EFG possesses axial symmetry. In this case, 
the quadrupole energy is completely defined by the spin 
projection I\ and the quadrupole quantization axis is di- 
rected along e\. 

For the concreteness we will particularize the case 
1=1 and the case of axial EFG symmetry x = 0. 
These assumptions simplify the problem noticeably. At 
the same time, it provides a good qualitative picture of 
the phenomena. In addition, some results will be pre- 
sented for the case I = 3/2. 

Let us choose the Cartesian reference system 
e x ,e y ,e z . Consider the eigenfunctions of the operator 
I x , \— 1) , |0) , |1) , as a basis. Then, the Zeeman inter- 
action, which is the same in both wells, reads 



H m = -gl3B\ ^ 







sin 9 , 
n/2 ' 



V2 



-UP 



V2 
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Here the angles <p and 8 assign the magnetic field direc- 
tion. 

Let e and e' be the axis of the EFG in the left and 
the right well respectively and a the angle between e and 
e'. Without loss of generality one can assume that e and 
e' belong to the x — y plane and = e. Then, in the 
right well 



Hr- 



3 „2ia 



|e" 2M 



(7) 
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while in the left well 



Hl = H Q (0) . 



(8) 



Thus, the Hamiltonian of the tunneling particle in the 
presence of the quadrupole and Zeeman splittings reads 



TT-l ( Hl Ao ' 1 
2 \ A • I H R 

where I is a unit matrix of rank 3 and 



Hl = -A ■ 1+2- (H LQ + H„ 
H R = A-I+2- (Hrq + H m .) 



(9) 



(10) 



Parameter b is supposed to be approximately indepen- 
dent on the position of the tunneling particle within a 
wells. Also, we assume this parameter is almost the same 
in different wells. This is justified by the experimen- 
tal data for nuclear quadrupole resonance in amorphous 
AS2S3 and As2Se3 where only about 10% line widening 
has been foundSi 

Almost nothing is known about the angle a between 
the quadrupole axes in the DWP except for the case of 
amorphous glycerol where this angle has been found to 
have a specific valued. The analytical approach below 
(see. Eq. I|31|0 reveals that there is an insignificant de- 
pendence of the result on this angle, perhaps excluding 
the cases where a — 0,7r/2,7r. For this reason, we as- 
sume the angle distribution to be uniform on the interval 
< a < 2?r. 

In this paper, for the sake of simplicity, we are re- 
stricting our consideration to the single tunneling par- 
ticle possessing a nuclear spin 1=1. This assumption 
simplifies the analysis. In general, higher spins and mul- 
tiple particles are possibly involved in a single tunnel- 
ing systemii2i2ii2£. We assert, however, that our consid- 
eration captures the important physics. In particular, 
our assertion is based on the fact that we have not re- 
vealed any qualitative difference between spins 1=1 
and / = 3/2. The effect of a multiparticle TLS struc- 
ture seems to be insignificant at high temperature where 
the effect of interest is a small perturbation and we can 
treat the nuclear quadrupole interaction for each particle 
belonging to a tunneling system independently. Then, 
the resulting effect is reduced to a sum of single parti- 
cle contributions. To take into account the complex TLS 
structure one should multiply the quadrupole constant 
& by a factor of the number of particles involved in the 
tunneling system. 



III. RESONANT PERMITTIVITY OF A MULTI 
- LEVEL TUNNELING SYSTEM 

In the model considered in this paper the application 
of the external electric field makes the asymmetry pa- 
rameter field dependent 



Here F is the strength of the field, and 



-I 
I 



(12) 



is the dipole-moment operator of the system described by 
the Hamiltonian ©, so that the dipole moment of the 
particle in the left well is — /^/2, while in the right well it 
is fi/2. 

In our recent paper— we have derived a general ex- 
pression for the dielectric permittivity for the ensemble 
of quantum system each having a discrete energy spec- 
trum Err, 



2n 



A'. 



ab = ^ ( -f+aHbP, 



d 2 E n 
' dA 2 



dP m dE m 
dF b dA 



(13) 



Here E m are the eigenvalues of the Hamiltonian @; fi a 
and Ft,— are the Cartesian coordinate of the vectors fi 
and F respectively; n = 21+1, e.g., n = 3 for the case 
of nuclear spin 1=1; P m is the population factor for the 
state E m . 

Generally, in the low frequency external electrical field 
F ~ exp (—2irivt) the population factor P m can change. 
This process, however, is important only when the tun- 
neling system experiences the relaxation for time r 
smaller then the period of electrical field oscillations v^ 1 . 
In the regime of interest, T < 50mK and v > 100-ffz,the 
time r is large enough in comparison with the field os- 
cillation period j/ _1 — Accordingly, the relaxation rate 
of TLS t -1 ~ v is so small that the population of en- 
ergy levels cannot follow the rapidly changing field F. 
For this reason, the population factor remains approxi- 
mately field independent— Therefore, the contribution 
of the second term to permittivity in Eci. (|13fl (usually 
called a relaxational contribution) vanishes. So the per- 
mittivity is only due to the first term in Eq. i|13|) called a 
resonant contribution. Then the population factor P m is 
given by the unperturbed (in the absence of the weak field 
F) equilibrium distribution for the by the field TLS's 



Pn 



exp (.— %*■ 



exp 



(14) 



The contribution of a single TLS should be summed over 
all TLS's belonging to the system. This is equivalent to 
averaging the susceptibility Xab in Eq. (|13f) over ener- 
gies, tunneling amplitudes, dipole moments, and angles 
between quadrupole axes and the external magnetic field 
(angles a,6,ip in Eqs. ©, 0). 

Taking the average over directions and absolute values 
of TLS dipole moments is straightforward and we can 
rewrite the resonant TLS contribution in the form 



2)1 



X r a,b = $ab--X, X = ( P ™Xmj , 

d 2 E m 



A (F) = A — F/x 



(11) 



X, 



dA 2 



(15) 
(16) 
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The parameter Xm is a function of tunneling parame- 
ters A, A , the external magnetic field B, and nuclear 
quadrupole interaction b in both wells. The angle bracket 
in Eq. I|f 5|l mean the Gibbs averaging, the averaging over 
distribution J2J), over directions of the magnetic field tp 
and 9 and of electrical field gradients a. Thus, the per- 
mittivity x remains the function of temperature T, the 
quadrupole interaction b and the magnetic field B. Note 
that one can ignore the external field F effect on the 
local electric field gradient interacting with the nuclear 
quadrupole because the field F is 5 to 6 orders of magni- 
tude smaller then the crystal field. 
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IV. PERMITTIVITY OF TUNNELING 
SYSTEMS WITH QUADRUPOLE 
INTERACTION. THE NUMERICAL ANALYSIS. 

Below we analyze the case of high temperature, i.e., 
the case when 

6<T. (17) 

It is known that the main contribution to permittivity of 
tunneling systems in glasses is due to tunneling systems 
having A ~ Ao > T4ii£ On the other hand, it is clearly 
explained iniS that merely the magnetic field (in the ab- 
sence of quadrupole interaction) does not change the per- 
mittivity. Changes can occur only if the quadrupole in- 
teraction exists. Below we argue that, in the absence 
of the magnetic field, the correction Sx to permittivity 
caused by the nuclear quadrupole interaction is associ- 
ated with tunneling systems having A < Aq < 6. 

The net quadrupole effect can be described by the dif- 
ference of two permittivities Sx, one given by Eq. I|15|l 
and the other defined by the same Eq. I|15|) taken at the 
value b = 

Sx = X-x(b = 0) (18) 

In fact, it is impossible to "turn off" the quadrupole inter- 
action to experimentally investigate how the quadrupole 
interaction influences the permittivity, if any. However, a 
very strong magnetic filed should inhibit the quadrupole 
effec1ji2i22i2i This is due to the fact that the strong mag- 
netic field creates identical local eigen states, charac- 
terized by a certain nuclear spin projection onto the 
magnetic field in the wells. Thus, the influence of the 
quadrupole effect on the permittivity can be verified ex- 
perimentally by measuring the difference 

= Xcx P {m = 0) - x(m = ao), m = gfiB, (19) 

where m is the Zeeman splitting of nuclear spin energy 
levels by the external magnetic field. 

The permittivity of the tunneling system is completely 
defined by the energy spectrum of the Hamiltonian 
This spectrum strongly depends on the relation between 
the nuclear quadrupole interaction b and the Zeeman 
splitting m. 



FIG. 2: A temperature dependence of the contribution to the 
permittivity of tunneling systems due the quadrupole interac- 
tion. The graphs clearly shows the 1/T dependence above the 
temperature T = 5mK. The graphs in the inset also shows a 
negative contribution below T = lmK. 



Consider first the case m <C b where one can neglect 
the magnetic field, i.e., approximately assuming m = 0. 
Let us first present the results of the numerical simula- 
tion for the permittivity (see. Fig. [21 EJ) in the case of 
b = lmif.These results were obtained as follows. First, 
we chosen a certain set of admissible values of the pa- 
rameters deterring the Hamiltonian @. Next, we found 
numerically the eigen values of the Hamiltonian and used 
Eqs. i|14|) . i|15|) . (|16fl to calculate the permittivity of a 
TLS with fixed parameters^ Then, numerical integra- 
tion of the single TLS responses over all relevant pa- 
rameters was performed to evaluate the effect of aver- 
aging. Fig. [21 represents the permittivity induced by the 
quadrupole interaction. For low temperatures T < 5 mK 
this correction is negative. Also, in this temperature 
range, as the temperature increased the permittivity also 
increased. These "low-temperature results" are due to 
the low temperature breakdown of coherent tunneling de- 
scribed in our previous papers^2*2i 

Above 5mK the quadrupole interaction induced per- 
mittivity falls with the temperature. The analysis reveals 
that this decreasing follows the dependence |6| /T. Figure 
121 displays this dependence neatly for both spins 1=1 
and 7 = 3/2. 

Using a similar numerical approach, one can investi- 
gate the magnetic filed dependence of the permittivity. 
In this case, the Zeeman splitting has been taken into 
account in Eq. Q with averaging the result over the 
direction of the magnetic field. Fig. [3] displays a sharp 
increase for the Zeeman splitting region m < 2b, then a 
quasi-plateau in the magnetic dependence takes place for 
2b < m < T followed by the decrease for 2 m > T. 

Let us now give the physical interpretation of these 
temperature dependencies. 
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FIG. 3: Quasi-plateau magnetic field dependence of the res- 
onant permittivity for the Zeeman splitting 5mK < m < 
20mK. The graph on the inset displays quasi exponential 
exp(— b/m) dependence for the Zeeman splitting m > 20mK. 



A. Analytical solution for the permittivity in the 

case 1 = 1 



b 



FIG. 4: The energy level structure for nuclear spin / = 1 in 
the vanishing magnetic field if the EFG possesses axial sim- 
metry. The thick solid lines denote the degenerate levels in 
the left and right wells. The thin solid lines denote the non- 
degenerated ones. The dotted lines connect the levels coupled 
by the tunneling abplitude Aq. 



There are four physically different cases m <§; b <C T, 
m<fe<T, 6<m<T, fe<T<m. 

First, let us consider the particular case, when the 
magnetic field vanishes. It was assumed and then con- 
firmed by the results obtained in this paper that only 
tunneling systems for which A, Ao < b significantly con- 
tribute to the correction 1(19(1 . For this case, we can 

assume E m <C T in Eq. ljl"4"|) and, therefore, e^"^ « 
1 — Since for the Hamiltonian © 



D 



0, 



(20) 



one can make the 1 [T expansion of the correction 

?7i (A,A ,6, a) 



T 

6 



m 



d 2 E v 
dA 2 



(21) 



Here the angle bracket mean the same averaging as in 
Eq. Ijl 5(1 excluding Gibbs averaging. 

For the Hamiltonian (jSJ one can find exact expressions 
for E m (see Appendix^). Also, the exact analytical aver- 
aging procedure is described in the Appendix. The final 
expression for the correction to the permittivity reads 



T 48 



0.2 • P'- 



T 



(22) 



This dependence shows a nonanalitic b- dependence 
which requires interpretation. 

Correction to the permittivity 8\ (b, T) induced by the 
quadrupole interaction b is due to tunneling systems hav- 
ing A Ao « b. It is evident that 5\ {b 7 T) can be es- 
timated as Sx (b — 0, T) taken for the tunneling systems 



having A « Ao ~ b. Thus the correction Sx (b, T) is esti- 
mated to be 

S X (b,T) 

)M 2 [ lb] dA fW Aid A 



P- 



3 Jo A J_\ b \ (A2 + A 2)3/2 



x tanh 



2T 



Pfi 2 \b\ 
T ' 



(23) 



This result correlates with the numerical analysis of Sec. 

El 



B. Qualitative estimate of S\ for intermediate, 
m rs \b\ , and strong, m> magnetic fields. 

Let us return to the case when tunneling systems are 
described by the Hamiltonian @ in the zero magnetic 
field. If one neglects the parameter A , the energy spec- 
trum in each of the well consists of one non-degenerated 
level and one double-degenerated level (see, Ref^ and 
Fig.H 

The contribution to the permittivity originates from 
the left and right states for which energy detuning is 
small compared with Ao. The tunneling term mutually 
bounds either the non-degenerate levels or the degener- 
ate levels. In this case there exist three different values 
of A that result in resonance. The magnetic field elimi- 
nates degeneracy and we obtain three pairs of mutually 
coupled levels and there appear six different values of 
A resulting in resonance. Thus, the effective number of 
tunneling systems contributing to the resonance increases 
(see. Refs-i^) 
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its contribution to the resonant permittivity a; 



,13,14 



A - m 



+1 



S\ (m,b,T) 

= P ti r ^ 

3 L A 



x tanh 



((A-m) + (A Ao ) : 



3/2 



Y(A-m) + (AA ) S 
2T 



P(i 2 n\b\ 
~~3 T~ 



(25) 



FIG. 5: The energy structure of a tunneling system in the 
case of a strong magnetic field b <C m < T. The resonance is 
shown to occur between the state with nuclear spin projection 
I z = — 1 in the left well and the the state with nuclear spin 
projection I z — in the right well. The energy mismatch is 
|A — ml <C m. 



In the case of a half-integer nuclear spin, energy levels 
in a zero magnetic field are degenerate according to the 
Kramers theorem. The application of a magnetic field 
eliminates this degeneracy. As in the case of integer nu- 
clear spin, this results in the increase of the number of 
resonances. This circumstance clarifies why the appli- 
cation of the magnetic field increases the resonant per- 
mittivity. Note that at T < b only resonance between 
the lowest levels in the left and the right wells is ther- 
modynamically allowed. So the effect described above 
vanishes ,22*21 

Let us then turn to the case of j&j <§; to. In this case, 
one can approximately classify energy levels by nuclear 
spin projection I z . If b — 0, the coupling between the 
left and right wells is possible only for levels having the 
same nuclear projection. If b is finite, the overlap factor 
between the nuclear spin states in the left and the right 
wells for levels with different nuclear spin projections is 
proportional to \b\ /m and can be approximated by 



*S" « ^A 
m 



(24) 



This expression results in the following remarkable con- 
clusion. In the case under consideration, the correction 
to permittivity is independent of the magnetic field. This 
conclusion does not completely correspond to the numer- 
ical simulation (see. Fig. |2J where a smooth temperature 
decrease in Sx was found. Therefore the result given by 
Eq. needs to be amended. While writing Eq. 
we did not take into account the population of the level 
with the nuclear spin projection I z = in the right well. 
This level is separated by an energy gap of the order of 
to from the level I z = — 1 in the right well which is the 
ground state of the tunneling system (see Fig. |5j. For 
this reason, the population of the level I z — in the 
right well is approximately exp (—m/T) . Therefore, the 
result given by Eq. I|25[) should be multiplied by the fac- 
tor exp (—m/T) resulting in a slow decrease of 5x with 
increasing m in the region 5mK < m < 30mK. Com- 
paring Eq. it22*|) and Eq. one concludes that the 
correction to the permittivity induced by the magnetic 
field is several times larger than for no magnetic field 
correction. This agrees with the numerical simulation 
presented above in Fig. |3 where the magnetic field first 
causes a sharp increase in the permittivity followed by a 
slow decrease. It is clear that the contribution to permit- 
tivity given by Eq. Q25[) is due to the transition between 
levels having different spin projections. This means that 
the contribution to permittivity is associated with the 
tunneling accompanied by a change of the nuclear spin 
projection. This effect is discussed below as a new po- 
tential mechanism for the spin lattice relaxation. 

If the magnetic field is strong enough that to 3> T, only 
the lowest Zeeman levels are occupied what brings us and 
we come to the case of the standard tunneling model. In 
this case, the quadrupole interaction effect disappears in 
agreement with Fig. [3] 



CONCLUSION 



Resonance between levels with I z = — 1 in the right well 
and I z — in the left well can be treated as a separate 
tunneling system having the energy mismatch A — to and 
the tunneling amplitude Aq^ w ^Aq. One can estimate 



It is shown in this paper that at temperatures T > 
lOmK the correction to permittivity induced by the 
quadrupole splitting b behaves as \b\ /T. This dependence 
correlates with the experimental data ML This compari- 
son with experimental data can be used to estimate the 
constant \b\. On the other hand, one can extract from the 
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nuclear quadrupole resonance (NQR) experiment the fre- 
quency of the transition lo . If the tunneling particle car- 
ries several atoms, the effect of the number of the atoms is 
additive within the perturbation approach. Then, the ra- 
tio \b\ /hujQ approximately estimates the number of atoms 
per a tunneling system. 

Our results exhibit a good qualitative agreement with 
experiment: there is a sharp increase in permittivity as 
the magnetic filed is applied, followed by slow decrease af- 
ter the Zeeman splitting passes the value of quadrupole 
splitting (see Fig. |3J). In fact, in BK7 glass a simi- 
lar effect takes place with position of the maximum at 
the magnetic field of the order of Teslas£ This agrees 
with the results obtained in the present paper. On the 
other hand, in BaO — AI2O3 — S1O2 glass the effect is 
observed at magnetic field about three orders of magni- 
tude smaller The reason for this effect is not clear. 
This can be due to the presence of residual paramagnetic 
impurities interacting with tunneling systems since the 
electronic spin in the magnetic field of the order of a few 
milliTeslas. acquires the Zeeman energy comparable to 
the temperature. 

However, one can suggest an alternative interpretation 
of the effect in BaO — AI2O3 — SiC>2- In Ref.— we have 
shown that resonant pairs of tunneling systems experi- 
encing quadrupole splitting noticeably contribute to di- 
electric loss (imaginary part of the permittivity). It has 
been shown that the behavior of these pairs resembles 
the behavior of two-level tunneling systems. The effec- 
tive energy splitting for them is of the order of 10/iA'. 
This value corresponds to the Zeeman splitting induced 
by the magnetic field of the order of milliTeslas. Thus, 
the magnetic field of the order of milliTeslas influences 
the behavior of tunneling pairs affecting the imaginary 
part of x by means of change in TLS relaxation rate^ 
Because of the Kramers - Kronig relation, this circum- 
stance should be reflected in the real part of permittivity. 

According to the experimenlj 2 *^, the position of the 
maximum in the dependence of the dielectric constant 
8\ on the magnetic field is sensitive to the value of the 
external electric field used in the measurements of the 
permittivity. This suggests the nonlinear character of the 
effect. The interpretation of the nonlinear effect has been 
proposed in Ref.™ It was found that the maximum posi- 
tion is proportional to the electric field. The dependence 
on the electric field qualitatively resembles the behav- 
ior in Fig. |31 In the present work, we investigated only 
the linear response. So, to attain the linear regime, one 
should experiment at lower electric fields. The indication 
to the linear response is the peak position independence 
of the electric field(see Fig|3J). 

It is interesting that the correction to the permittivity 
at large magnetic fields m > b is caused by tunneling 
accompanied by a change in the nuclear spin projection, 
Fig. [5] We are justified to anticipate that at low tem- 
perature a similar processes can contribute to spin - lat- 
tice relaxation with the rate exceeding that for dielectric 
crystals by many orders of magnitude. This new spin 



- lattice mechanism should be sensitive to the resonant 
transition rate. We expect that this transition is induced 
by spectral diffusion 28 which makes the transition an ir- 
reversible one. This relaxation can be slow down by the 
nuclear spin diffusion needed to bring a nuclear spin to 
the TLS neighborhood^ The analysis of this effect is 
beyond of the scope of this paper. 
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APPENDIX 

For the Hamiltonian JjJJ one can find the exact expres- 
sion for E m . The eigen energies of the Hamiltonian © 
in the case / = 1 at zero magnetic field read 



(26) 

E = v /(A/&) 2 + (A /fe) 2 ,M = 1 + E 2 , 
Y = 2VN, N = E 2 - (A jbf sin 2 a. 
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Substituting these expressions into Eq. I|21[l after simpli- 
fication one obtains 



Vi (A, A ,6, a) 
(A/6) 2 (A /6) 2 (E 2 + 1 - 4JV) sin 2 a 
6£ 2 iV (47V - (S 2 + l) 2 ) ' 



(27) 



Note that, if a = or a = tt, the quadrupole effect 
disappears in agreement with Ea. i|27[) . 

Then the total permittivity is given by the expression 



T 3 2 



it A J-T/|fc| 

x?7i (A, Ao, b, a) . 



(28) 



To calculate the last integral over A, Ao it is convenient 
to change the variables as follows. 



A = E ■ sine, A = E ■ cose; |e| < tt/2 
Then the correction H28JI takes the form 

rn {E, e b,a) = 6g2 , (4 ^_ M2) 
g = V 1 — c 2 , c = sin a ■ cos e. 



The indefinite integral over a parameter E reads 



(29) 



(30) 



r)i (E, e, 6, a) cLE 



— sin e • arctan 
12 ^ 



2- E-c 
1-E 2 



To calculate the integral (|28|l . one can substitute the 
integral limit T j |&| by oo because of the fast integral con- 
vergence. The first term in the rhs of Eq. 113 111 is a broken 
function of the argument E at the point E = 1. There- 
fore, integration should be performed independently over 
the two intervals (0, 1) and (1, oo) . After integrating over 
E one obtains the intermediate result 



Vi (e, ") 



12 



sin e • 1 



[ ^ ^ ' — - de — — | cot a ■ ( arctan ( 
J cose 12 V V V 



2-c , [2-E-g 

• arctanh — 

9 \1 + E 2 



r 



The contribution of the second term in the rhs of Eq. 
(|31|l into the integral (|28l) vanishes since it turns to zero 
on the both limits E = 0, oo. 



The indefinite integral over e reads 



tan e \ / sin e 

2 arctan 

cos a 



1 + cos £ 



cos a — sin e • cos e • sin a 



and 



Vi (") 



12 sin a 



cose 



rfe 



(1 — | cosa| 



Finally, after the integration over the angle a one obtains 



Sx (T,b) lx 



(31) 



T 48 



(4-tt) w 0.2 
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